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Abstract 



We calculate dispersion relations and propagators for the cr, w, tt, p, 6, rj and ai mesons in 

relativistic, dense, hot, symmetric nuclear matter. In addition to the usual mixing of the a — uj 

system, we obtain mixing of the S with the longitudinal p mode and of the tranverse p with 

the transverse ai mode. Finally, the component of the oi polarization along the transferred 

momentum modifies the in-medium pion propagator in a way similar to the Migdal contact 

^ ' interaction, but with the opposite sign. The spurious pion condensate as well as the additional 

'^ . contribution from the oi meson are removed by a contact term. We compare two ways of 

3^ '_ implementing contact term subtraction. 

,_i ; 1 Introduction 

^^ , This works investigates some points concerning the dispersion relations and effective masses of 
mesons in dense, hot, symmetric nuclear matter. The dispersion relations will be calculated in the 
^ ' random phase approximation (RPA), by performing a linear response analysis 0]-p| around the 
1^ . homogeneous Hartree ground state. This paper reproduces previous results and extends them in 
two aspects. 

Equations are presented for the six mesons of the Bonn model ^, i.e. a, to, it, p, 6 and ry. In 
^ . addition, the contribution of NN loops to the dispersion relation of the ai is investigated, since this 
H I meson can be of interest from various points of view. One is the construction of chiral Lagrangians, 
since the oi is the chiral partner of the p meson [^ - [O]. Another is the issue of the dilepton 



production in relativistic heavy ion collisions. As a matter of fact, the ai [16| couples to the p 
meson and may modify its spectral properties. The ai may also be used to describe part of the 
nucleon-nucleon interaction which is mediated by correlated n-p exchange in the S-wave channel 
[17, |l^. The inclusion of correlated n-p exchange is known to improve the Bonn potential p^ ; 
however one should keep in mind that the structure of correlated n-p exchange is more complex 
than what can be described by exchange of particles with a sharp mass ||l^. Finally, the ai meson 
has sometimes been quoted ]2y, ^] to provide a means to improve the behaviour of the differential 
nucleon-nucleon cross section for backwards scattering angle. This is related to the behavior of the 
tensor part of the NN potential near r = 0. 

As a first part of this work, a discussion of the ai meson dispersion relation and mixing effects 
with the p and vr mesons is presented in section §2. It will be seen that the transverse mode of the 
oi mixes with that of the p. Since the axial current is not conserved, we will have a part of the ai 
polarization proportional to q^q'^, where q^ is the transferred quadrimomentum. This part mixes 



with the pion and renormahzes its polarization in a way which is formally similar to that produced 
by a contact interaction of the Landau-Migdal type. However, it has the opposite sign, so that the 
oi would enhance the unobserved pion zero sound mode. Since experimental evidence shows that 
pion condensation does not occur at saturation density in symmetric matter, this spurious mode 
should be removed. The standard procedure to implement this is to introduce a contact interaction. 

We are thus lead in a second part to make some considerations on the relativistic generalization 
of the contact term. We will compare in section §3 the results obtained using the Ansatz of 
Horowitz et al. p2| , [2^ , who modify the pion propagator in an ad-hoc way, to a more standard 
implementation simply consisting of adding one or more contact pieces to the Lagrangian density. 
It will be seen that both procedures may lead to different results, when the former is not used 
with some caution. In particular, this may explain the discrepancy recently obtained in calculating 
the RPA corrections to neutrino- nucleon scattering in proto neutron star matter [^, ^, [2^. We 
advocate for the latter procedure. 

Finally, numerical results are presented in section §4. In symmetric matter, the dispersion 
relation factors in several subsystems. First there is the a-iv sector, with a transverse oj mode 
and a longitudinal mixed a-uj mode. Since this mode has already been studied extensively in the 
litterature, no numerical results are shown here. Then we have a d-p-ir-al sector, which is composed 
of three parts: a longitudinal mixed 5-p mode, a transverse mixed p-ai mode and a mixed vr-oi 
mode. Dispersion relations are calculated for these modes and the behaviour of the effective masses 
of the mesons is investigated. We mention the problem of vacuum fluctuations and of the choice 
of a renormalization procedure. As already noticed in the case of the p meson Q, different results 
are obtained depending of the choice of the renormalization conditions. At the time of writing, we 
must regard this as an unresolved issue. 

Among the possible applications of the present results, one can quote the calculation of screened 



nucleon- nucleon potentials [||, ^] or cross sections |28, 29| and their application to nuclear dynamics 



Ol, dilepton production in relativistic heavy ion collisions [10], the calculation of the electronic re- 



sponse function [^] or of the compressional modes of nuclei [^, 33, 34|. Another field of application 



concerns the physics of supernovae explosions and protoneuton star cooling. In conditions of high 
density and temperature reached in the early phases of neutron star formation, it has recently been 
argued that RPA corrections to the neutrino-nucleon scattering cross section could be responsible 



for a sizeable modification of the neutrino mean free path |24, 25, 26]. The results obtained in this 
work may be applied in the limiting case of pure neutron matter, by appropriately modifying the 
degeneracy factor. In astrophysical applications nevertheless, one has to deal with an asymmetric 
environment. The dispersion relation of neutral mesons in asymmetric nuclear matter is presented 
in a companion paper ]p5[. The case of propagation of charged mesons in asymmetric matter is 
currently under consideration. 

2 Meson exchange model 

2.1 Discussion about the mesons included in this work 

The a, uj, p and vr mesons are standard mesons which are included in the description of the NN in- 
teraction in hadronic models, whose prototype is the Quantum Hadrodynamics (QHD) developped 



by Walecka and coworkers ]36]. A vast amount of literature has been devoted to study their dis- 



persion relations ]U~lg], ]37]-]44]. This section presents a discussion on the relevance of extending 
the approach to study more mesons. 

The 6 and rj mesons belong to the set of mesons exchanged in the Bonn potential JH]. However, 
they have been less studied than the previous ones. One reason is that the authors of ]p found that 
they only bring a small adjustment to the form of the NN potential. In mean field studies of nuclear 



matter, the 6 is usually not considered on the argument that the exchange of p mesons is enough 
to reproduce the asymmetry energy at the mean field level (see however [45, ff^])- In asymmetric 



matter, on the other hand, the delta meson could be important, since it carries isospin. Interest has 
been shown for the S meson in the context of Dirac-Brueckner calculations of the equation of state 
in asymmetric matter. When developping equivalent mean field theories with density dependent 
couplings, in order to reproduce the results of the full many-body calculation, de Jong and Lenske 



[46 1 or Shen et al. [^] have found that it was necessary to introduce the 5 meson with a significant 
coupling strength, e.g. 5f|(nsat)/(47r) = 4.61 at saturation density. 

The S meson mixes with the p meson, just as does the a with the uj in the isoscalar sector. It is 
well known that a-u mixing is very strong (see e.e [^, ^), so that one may ask whether the same 
occurs in the case of 6-p mixing. In relation with the debated case of the interpretation of dilepton 
measurements in heavy ion collisions, the effect of 5-p mixing has been studied by Teoredescu et 



al. [44|. 



The oi is very massive (ma=1260 MeV), and is for this reason usually not included in meson 
exchange models, which take 1 GeV as a reasonable cutoff energy scale. There are however theo- 
retical reasons to consider this meson. As mentioned in the introduction, there is a need of keeping 
this meson from symmetry arguments since the oi is the chiral partner of the /O [0] - p!5| ]. 

For the purpose of reproducing the dilepton measurements of CERES and HELIOS, an enormous 
amount of work has been spent on calculating the p spectral function, as this meson yields the 
dominant contribution to dilepton production in the vector dominance model. The controversy 
arised, as to whether the dilepton enhancement at lower invariant mass was due to a reduction 
of the p meson mass in the medium and a signal of chiral restoration in the Brown- Rho scaling 
picture [|4^ (B/R scenario), or could it be explained by purely hadronic models as an increased 



width and shifted strength due pmr coupling [49[ (R/W scenario). It was suggested [IT^, ISOl that 



one may expect a substantial contribution to the width from pnai loops. 

An other ground to study the oi meson is to determine whether it could simulate short range 



corrections in the pion exchange potential [^, 21 1. As a matter of fact, Fourier transforming the 



expression obtained for the pion exchange contribution to the transition matrix yields a 5{r) piece 
in the NN potential in coordinate space. The standard procedure is to introduce the so-called 
Landau-Migdal parameter g' as a contact interaction to remove this piece. This term is necessary 
to achieve a good description of the proton-neutron cross section at large scattering angles 9 ~ 180" 



[51, ^]. The parameter g' then also enters in the expression of the pion polarization and further 



has the nice property of removing spurious zero-sound modes from the pion dispersion relation. 



which would be responsible for the onset of pion condensation at unrealistically low densities [53]. 

In a relativistic formalism, it is less clear how to implement this term in a simple fashion. 

For example, Horowitz et al. [^, ^ have suggested an Ansatz making a replacement in the 

pion propagator. An other possibility would be to add to the Lagrangian a contact term with 



pseudovector coupling. This is the method used by Schafer el al. [54[, who add to their Lagrangian 
a term gai'4^75l^T"ip)itp75^fj,T"4') in order to fit the NN cross section at large scattering angles. 
Equivalently Engel et al. [^] suggest to add to the Lagrangian a piece gaV'757'^^o^ and take the 
limit of infinitely heavy a^ field. In section ^, we compare the effects of the Ansatze of Horowitz et 
al. or of Schafer et al on the one hand, and of true oi exchange on the other hand, on the dispersion 
relation of the vr and p mesons. It will be seen that the mixing of the oi meson with the pion acts 
as a Landau-Migdal term with the "wrong" sign. 

2.2 Outline of derivation 

We will take the following Lagrangian density 

C = /^ATTV + ^Af$ + ^<J>$ + >C(j)j$2 + /2cT (1) 
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The first two hnes are the Lagrangian for the nucleon and its interaction with the meson fields. 
We chose the pseudovector coupling for the pion since its phenomenology is better, namely it 
permits to reproduce the pion-nucleon scattering. The pseudoscalar coupling will also be considered 
in section 



4.4 



The following two lines contain the kinetic terms for the mesons. The last two lines contain 
some meson-meson interactions which are important from the point of view of phenomenology. In 
particular, we may have frvrvr, pyrvr and pvrai couplings, which contribute to the spectral width of 
the o", p and ai mesons. The sigma self couplings bma"^ + ca^/^ improve the description of the 
effective nucleon mass and incompressibility modulus. 

One could consider more complicated forms of the Lagrangian, which would be obtained by 
enforcing chiral symmetry pi O], [l^ - |14|. Accordingly, more meson- meson coupling terms 
would appear in £$^$2- Since our purpose here was mainly to investigate NN loops, we will not 
consider such terms in this work. In fact, as will be seen in the following, the derivation method 
based on linear response analysis, which we follow here, discards meson-meson correlations at an 
early stage of the calculation, so that terms involving meson couplings only appear at the mean field 
level. Since we deal in this paper with symmetric matter, most of these terms actually disappear, 
since they would involve expectation values of the pion, rho or ai fields, all of which vanish in this 
approximation. Only terms involving the a field will survive. Meson loops should of course appear 
in a more realistic treatment. 

The term Cct represents a counterterm Lagrangian in order to handle the vacuum divergences; 
its form is given in Appendix A. 3. 

The dispersion relations are derived from a linear response analysis in the Wigner operator 
formalism. By applying a perturbation to the relativistic Hartree equilibrium, one obtains homo- 
geneous equations for the perturbed meson fields 4)i{q) in the form D{q)(j)i{q) = 0, where (l)i{q) is a 
vector formed by the components of the perturbations to the various fields a, uj, it, p, 6, rj and oi, 
and D{q) is a matrix which contains the dispersion relations. The condition that these homoge- 
neous equations admit non trivial solutions is that the determinant of D(q) vanishes. As D{q) is in 
general non diagonal, the determinant will only partially factorize. A standard example is the case 
of the a-uj subsystem, where the uj transverse mode factorizes out, while the longitudinal mode of 
the u! mixes with the a mode |2| It was shown in this work that the results obtained coincide with 
the one-loop approximation from Green's function formalism. 

We define the Wigner operator 
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It obeys the kinetic equation (as well as a conjugate equation) 






with $(x) defined as in Eq. @ At this level, F{x,p) as well as the fields o-{x), uj{x), tt{x), 6{x) ...etc 
contained in ^ are operators. In order to obtain the RPA approximation, two basic assumptions 
are made: 

— Correlations are neglected. When taking the statistical average of Eq. (^, it amounts to 
replacing in the right hand side the average of the product < ^F > by the product of averages 
< ^ X F >. It is at this point that we neglect the contribution of meson loops to the dispersion 
relations. They could be restored by releasing this assumption, or be reintroduced by hand at the 
end of the calculation. We will not consider them in this work. We will omit the <> denoting 
statistical averages henceforward in order to simplify the notations. 

— It is assumed that there exists an uniform, unpolarized equilibrium given by the relativistic 
Hartree approximation, and that one may perform an expansion around this equilibrium (remember 
the notations F, a ... represent from now on the statistical averages < F >, < a > ...) 

F{x,p) = Fh{p) + Fi{x,p) , a{x) = an + cri{x) , 5{x) = 5h + Si{x) 
a;^(x)=< + u;f(x) , f^{x) = ffj, + ff,{x) (7) 

tt{x) = TTi{x) , 7]{x) = r]i{x) , a'^{x) = Si^x) 

In eqs. (|7|) the pion, eta and ai contributions vanish due to parity arguments in unpolarized matter. 
Moreover, in symmetric matter the chemical potentials of the proton and neutron coincide, so that 
Ph and Sh will vanish as well, since they are proportional to the difference between the proton and 
neutron density, and between the proton and neutron effective masses respectively. (We do not 
take into account the tiny difference due to the electromagnetic interaction.) 

With these approximations, we obtain after linearizing and Fourier transforming Eq. (|6|) the 
first order perturbation to the nucleon Wigner function 



with 



Fiiq,p) = G[p-^] <l>{q)FH[p+^]+FH[p-^] <^{q)G[p+^] (8) 
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+9pl^Pi^.{q)-r + ^(T^''q^PiAq)-r + gaj5l^aif,{q).T (9) 

7.P + M 

Fh{p) = S{p)ip{p) 

with S{p)=-f.P + M (11) 

and ^{p) = j^ 6iP^ - m2) [e{po)n{p) + d{-p^)n{p) - 9{-p^)] 

In symmetric nuclear matter, d = 2 is the isospin degeneracy. n[p) and n{p) are the Fermi-Dirac 
distribution functions for the (quasi)-particles and antiparticles respectively. M = m — gaCn is the 
effective mass and P^ = p^ — g^jU)^ is the effective momentum. Finally, we insert this solution in 
the linearized equations of the mesons. For example, we have for the ai meson 

qi^q- + (_g2 + ^2)^M-] ^^^(^) ^ -gap^^i{q) X ff^ + ga J d^ pTr[^5l''rFi{q,p)] (12) 



Since we assume in this work that we are in symmetric matter, the Hartree component of the p 
field vanishes. Inserting ^ in this expression, we can recast it in the form 



-q^q"" + {q^-ml)g>^-\ ai.{q) = K^{q) Mq) +Kp{q) Piu{q)+Ka{q) aiu{q) 
where, for example 
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Terms such as Hauj vanish in symmetric matter since the trace over isospin will involve the difference 
of proton and neutron distribution functions. The dispersion relations may be summarized in 
matricial form D{q)4>{q) = In symmetric matter the dispersion relation decouples in three blocks: 
(cj-u;), (r/) and {p-d-ai-Tr): 



with 
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(15) 



In the preceding equations, the meson self couplings were not taken into account. If they are 
present, they appear as mean field modifications. The dispersion relations in this case are given by 
making the following replacements for the a and vr masses 



m„ 



Ml = ml + 2bmaH + Sccth 
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We will decompose the polarizations on the usual orthogonal set of tensors and vectors as follows 
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The explicit expression of the polarizations can be found in the Appendix. 

We obtain the dispersion relations by equating the determinant of the matrix D to zero 
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2.3 Propagators for the a-uj-Ti-p-d-'q-ai model 

The propagators G are obtained by inverting the dispersion matrix D. 

The a — oj sector has already been studied thoroughfully (see e.g. [||, |^, ^, |3|, ^, ^, Q). We 
obtain as usual 
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The a meson mixes with the longitudinal part of the uj meson. The mixing effect is known to be 
strong. A zero sound mode may appear in the longitudinal mode, with a strength which depends 
on the values of the cutoff parameters in the nucleon form factors [^, ^^ . 

The r] meson decouples from the other ones in symmetric matter, therefore the calculation of 
its propagator is trivial. We have: 
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If we had not taken into account the ai meson, the p-6 system would decouple and reproduce the 
pattern of the a-uj in the isovector sector. Now the p also couples to the oi through the polarization 
li-ap- Nevertheless, the calculation shows that the dispersion relation of the (p-S-ai-ir) subsystem 
factorizes in three terms: a mode where the 5 meson mixes with the longitudinal part of the p, a 
mode where the transverse contribution from the ai mixes with the transverse part of the p, and 
a mode where the pion dispersion relation is modified by a term coming from the part of the ai 
meson which is parallel to q^q'^ . We have 
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The mixing with the oi meson acts as a modification to the pion polarization. Using the 
relations existing between 11^, Ha-n and Hag (see Appendix A), it can be recast into the form 
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It therefore comes out in a form similar as would a Landau-Migdal interaction. We would obtain 
exactly the Landau-Migdal form in the relativistic case with the replacement 
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We must note however the sign in the previous equation. One would need a fictitious ai field with 
an imaginary coupling so that equation (|40| ) may be fulfilled. Contrary to some expectations raised 
in the literature |2^, the oi meson cannot be used to improve the short range behaviour of the 
pion-nucleon interaction; instead it contributes an additional term which must be compensated by 
e.g. a residual contact interaction. 

Also when calculating the potential generated by the exchange of vr and ai mesons in vacuum, 
one obtains after taking the semiclassical limit k/m <C 1 and also k/rria ^ 1 
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and the term responsible for the 5{r) singularity (Ti.(T2 after Fourier transforming to position space 
would be removed by the same choice [g'^/m'^) —>■ — (1/3) (/^/m^). 



3 Short range behavior 



As is well known, and also will be seen in section § [4.4| , the dispersion relations derived in the 
preceding paragraph would lead to an excessive softening of the pion mode and to pion condensation, 
which is not observed experimentally. This is due to the fact that we have not yet taken into account 
the effect of short range corrections at this level of approximation. A related shortcoming of the pion 
exchange model is that it gives a vanishing differential cross section at scattering angle 9 = 180° in 
the neutron-proton exchange reaction, in glaring contradiction with the experiment. Both features 
are due to the fact that the contribution to the NN potential arising from pion exchange Vn{r) 
contain a piece (5(r) singular at the origin. This is an artefact of the model due to the assumption 
that the particles are pointlike. Such singular pieces in fact appear for all kinds of meson exchange. 
Whereas it is enough for the a, oj terms entering the definition of the central potential to smooth 
this divergence by folding it with a form factor of the type 




a e {o-,a;,7r,p, ...}, 
it is necessary in the case of the pion to remove the 5{r) function in order to arrive at a correct 



description of the NN cross cross section and spin transfer observables [51|. 



Short range corrections of the desired type would arise in a many body calculation from higher 



order correlations at the ttNN vertex |56|. In a simpler approach, the usual practice is to Fourier 
transform the potential to coordinate space and there remove the delta function by hand. This is 
equivalent to adding to the potential the Migdal contact term g'cri.a2Ti.T2. The delta substraction 
procedure for the pion alone amounts to taking g' = 0.333. A further contribution comes from the 
rho. Phenomenology favors higher values g' = 0.5 — 0.9. A recent analysis [^ extracted g' = 0.6 
from data on the Gamow Teller resonance. This procedure works well in the non-relativistic limit 
[58 1 but has the disadvantage of not beeing covariant. 



In the following we examine two procedures which have been suggested in the litterature to 
introduce the Landau-Migdal short range correction in a covariant way. 

3.1 Propagator replacement Ansatz 

Horowitz et al |^, ^] suggest to replace the vertex and propagator of the pion as follows: 
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In the medium, the pi and rho propagators obey the coupled Dyson equation G = G^ + G^HG or 
equivalently, [G]~^ = [Go]~^ — 11, or explicitely 
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The "pion" polarization which appears in this expression is given by 
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and the usual pseudovector pion polarization is recovered by contracting with q^q^ 

^pv ^ _nM-g^g^ = -q^ii^Q (46) 



HZ = -(^) /dVTr 



ihi^s{p - -)-/5YS{p + -) 



(44) 



This "pion" mixes with the p meson, due to the polarization 



m". 



TTp 



id pTr 






757^^ ) S{p - I) [ gpY + {^'y''\x ] Sip + ^) 



ip{p - ^) - ip{p + ^] 



2p.q — ie 



QpVx 



Det(G-i^ 



m. 



(^ + n,T)(<?' -ml + UpT) - r?'g'n^. 



n2 



2 2 

q -m^ 



By inversion of (|43|), we obtain the propagators 



X (l + (7'n^L) X - 



c/'[g2 - g'{ml - q^)] 



Q{H)t,U 



-G[fL>'^ - G^lpT'"' - cf^Q^''' 



G 



'pL 

{H) 
pL 



1 



G 



q"^ -mj + UpL 
(H) L 



n2 

^^pTT 



q^ -m? + lipT - vi^q^g' , — 



^-' PTT 



(^{H)f^u 



^PQ ^2 

Gp^ e'^'P^pm 
G(^) = 



-n 



pTT 



' pTT 



2„2tt2 



(g2 _ ^2 + npT)(- + n^T) - v'q'^ 



pTT 



.(^) 



.(^) 



(H), 



f-iK"^) jpv _ ^(^)rppu _ ny'^lQ. 



)PV 



G. 



(H) 
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1 + g'n^L 



G 



(H) 



(g2 -ml + UpT) 



2 2tt2 



G 



ttT ~ 1 

(g2 _ ^2 + n,T)(^ + n^T) - ^'g'n^. 

(H) [g'jml - q^) + g^] 



{ml-q^){l+g'Ii^Q)+q^Ii^Q 



(47) 



The full dispersion relation in the vr-p sector is obtained by taking the determinant of (|43|). We 
have 



(48) 



(49) 
(50) 

(51) 
(52) 



(53) 
(54) 

(55) 



Using the relations between the polarizations Hap, ^-np, ^aL, n^^, it can be checked that one 
can go from the equations for the p-ai-ir system to those resulting from the Ansatz of Horowitz by 
taking the limit q^ /m^ -^ and replacing {ga/'nia)'^{m-K/ f-n)^ by —g\ so that we see again that the 
mixing with the ai meson acts as a Landau-Migdal term with the opposite sign. 



3.2 Contact term 

The procedure consisting of adding a contact term to the Lagrangian is straightforward 



C3i) 



- d -m —^5l'^dnTT.f 

2 m^ 



V' - 9Ai'>p75ltJ.Ti))-{i^l5YTi:) 



(56) 
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The new term does not directly modify the field equations for the pi meson. Rather, it contributes 
an additional pseudovector term to the evolution equation for the Wigner function. Following the 
derivation method used above, after neglecting all correlations, we obtain at first order 



T^PM-f)-"^ 


Fi{x,p) B 


-q^ + ml 


TTi = i 





fn 



rrij 



l5l^qiMr.7:i{q) + 2gA-f5l^T.Pi{q) 



with 



Pf(g)= d^kTT[^,rrF,iq,k)] 



Fh{p+1) (57) 
(58) 



(59) 



Taking the trace of Eq. (^), where ^{q) should now be replaced by ^{q) + 2(7^757^fPC(g), with 
757^^ and integrating, yields a self consistent equation for P^ 



priQ) 



i-^WZq^T^iiq) + 2gAllAAPiu{q) 



(60) 



The polarization 



n 



l5rrGip - |)757'rF^(P + |) + (G ^ Fh) 



= / d'^p Tr 
has already been met before. The equation is easily solved as 



Piqv = T'\ — n- 



TTi 



(61) 



Replacing this solution in the first order term of the field equation for the pion, we obtain the 
dispersion relation of the pion modified by a contact term 



m. 



U 



m-n 



q^^\ 



^J l-2gATl 



Q 



(62) 



With — (/7r/"i7r)^g^nx4 = IIpv, and defining 2gA{'m.,r/fTr)'^ = g', we recover the standard expression 
for the pion pseudovector polarization modified by the Landau-Migdal g' parameter. 



n. 



f. 



m-n 



q'm 



Hp 



1 - 2gATriA r - ff'Hp 



(63) 



Apart from redefining the pion polarization, the gA contact term also mixes with the transverse 
part of the p meson. The result is given in Eq. (|73[). Let us nevertheless study before a more 
general choice of the contact interaction. The rho meson exchange potential also gives rise to a 
delta function in the coordinate representation of the potential (see eq. (|114| )). This delta can be 
eliminated independently from the one appearing in the pion potential if we add another contact 
term for the p meson 



C3i) 



F 



m - gpTPf,- 



-R^,R^' + -mlp^p^ - 5i?(^7M^^)-(V'7''rV') 



(64) 
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When solving the equation for the Wigner function, there is an additional contribution to the 
isovector density. At first order of the linearization procedure, we have 



Q 



Fi{x,p) 



with 



gpl'^Tpi^.iq) + 2gR-f''T.Rif,{q)\ Fh{p + |) + 



R'liq)= d^kTr[rTFiiq,k)] 



(65) 



(66) 



In symmetric matter, the mean value of the rho Hartree field and the contribution to the isospin 
J d'^p Tr[y>^TFH] vanish. Ri{q) obeys the self-consistent equation 



current /^ 



Ri{q) = -gp^vvPiv 



29rTI^^R 



Iv 



(67) 



where Ilyy is the vector contribution to the p polarization IIjJ^'' = ^pllyy + 2gpfpI[YT + fp^TT- It 
is solved by projecting this relation on T^'^ , L^^^ and Q'^^ . In the general case that pion, rho 
and oi fields are present, and we moreover introduce both the contact terms 5yi('0757/i''"^)^ ^'^id 
gRii^l^T'4')'^-: the various conponents obey coupled equations. As before, the longitudinal modes 
decouple, whereas the transverse mode of the oi is mixed with the p and the pion mode is modified 
by the part of the a\ polarization which is parallel to q^^q'^ . We arrive at the dispersion relations 



m-™ 



u 



n2.g2 



(g2 _ ^2 ^ npr)r^Vi. + Hp^TT^^ai, = 
{q^-ml + fipLWpi, = {) 



vfi(<?) = 



(68) 

(69) 
(70) 
(71) 
(72) 



The explicit form of the polarizations corrected by the contact terms qa, Qr is given in Appendix 
B. When gR=0, the polarizations of the p reduce to 



HpT 
flpL 



IlpT + 



2gA{IlAu)^ 

i-25AnL 



n 



pL 



(73) 



The equations (|63|,|73D are in fact the same as obtained from Horowitz propagator replacement 
Ansatz (see the third factor in eq. ( |48| ) and eq. (|50D). One should notice however that the later 
Ansatz, apart from producing the desired result for the pi and transverse rho propagators, also 



introduces additional spurious terms in the dispersion relation (48) arising from the unphysical 
part of the new "pion" propagator orthogonal to q^q'^ like e.g. 1 + g'Tl-KT and 1 + g'H-KL- 



4 Dispersion relations: numerical results 

In this section we study the dispersion relations which result from the models described in the 
previous sections. We will concentrate on the S-p-ir-ai sector, since the a-LO branches have been 
studied extensively elsewhere. 
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4.1 Model parameters 



When not explicitely stated otherwise, the parameters are chosen to be those of the Bonn potential 
[0]. They are summarized in the fohowing table. 



meson a 


TJla [MeV] 


•2 
9a 

47r 


A„ 


a 


550. 


8.2797 


2000. 


UJ 


782.6 


20 


1500. 


vr 


138. 


14.6 


1300. 


P 


769. 


0.81 
( fp/gp=G.l) 


2000. 

{rip =2) 


5 


983. 


1.1075 


2000. 


V 


548.8 


5 


1500. 



The coupling of the pion is taken in this Bonn model to be g'^/{4:Tr) = 14.6. One usually prefers 
now a somewhat lower value 5f^/(47r) ~ 13.67 so that /t^ = mT^gT^/{2m) ~ 0.965. An other issue is 
the value of the cutoff, which is taken to be 1300 MeV in the Bonn model, whereas most calculations 
and experimental determinations of this parameter would give K^^ ~ 800 MeV. 

Data for the coupling of the ai meson to the nucleon is scarce. A simple realization of the 
chiral symmetry is Qa = Qp- From early chiral models by Weinberg and Wess & Zumino |0] one 
obtains the relation ga = 'ma{f-K/'m-,^). With m-a=1260 MeV and /tt — 0.965 as determined from 



NN scattering data, this yields ga — 8.79. More recent implementations of the chiral model |11, [I^ 
point out the difficulty in adjusting all known data on the width and mass of the ai and quote 
values ranging from 3.8 to 18 [15|. Another value used in the literature is ga = 6.44 [^]. Reference 
[^ also provides us with a rather low value for the cutoff AaNN = 809 MeV. Here we will work 
with ga = 8.8 and a higher value of the cutoff AaNN = 2000 MeV. 



4.2 Longitudinal d-p mode 

The dispersion relations were calculated with the parameters of Machleidt's Bonn potential Ig using 
three alternative renormalization schemes. As discussed e.g. in |5|, where the dispersion relation 
of the rho was studied without mixing, the polarizations contain a diverging contribution from 
vacuum fluctuations which has to be renormalized. One may define several subtraction schemes, 
which, unfortunately, lead to very different behaviors of the mesons effective masses. In ||5| two 
classes of renormalization procedures were identified 

^ a first class which renormalizes divergences of the form M^/e by subtracting a general coun- 
terterm Aq + Aia + ^2C^, and adjust the constants so that the polarization 11 and its derivatives 
(dli/da), {d'^n/da'^) vanish at some point. In this way, one tries to minimize the effect of introduc- 
ing new couplings with the sigma field in the counterterm Lagrangian, which are not present in the 
original physically motivated Lagrangian. In references |g, |3|, the renormalization point is chosen 
to be the mass shell q"^ = m^ for the polarization as well as its derivatives. In ||3^ Kurasawa and 
Suzuki subtract the polarization at q'^ = rn^ but the derivatives at q^ = 0. Using this procedure, 
the Lo meson mass decreases with density, however the p mass increases due to the tensor coupling 

/p- 

— a second class which preserves the structure of contributions Af ^/e = {m— g^a)'^ / e by subtracting 

a countertern in the form A{m—gcr(y)'^. A is determined by setting yl = at some point (at q^ = m^) 

Since we have only one parameter A, it is not possible to minimize the new couplings to the a field 

introduced in the counterterm Lagrangian. Using this procedure, the p meson mass decreases with 

density. 
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The reader is referred to ||5[ and Appendix A. 2 for further details. Once a renormaUzation 
procedure is chosen, it should be used for all mesons. 

The first renormalization scheme (scheme A) used in this paper belongs to the first class for all 
mesons. For the p meson, it is the "scheme 2" described in the Appendix C of P]. For the 5 meson, 
it is obtained by replacing everywhere m(j and (70- by ms and gs in the expression for the vacuum 
polarization of the a given in Ref. |^ 

The second renormalization scheme (scheme B) belongs to the second class for all mesons. For 
the p meson, it is the "scheme 3" described in the Appendix C of |^. 

A third renormalization scheme (scheme C) is used for comparison with calculations published 



in the litterature |44|. For the rho meson, it uses as in |44| the procedure of Shiomi and Hatsuda 
[42 1, which according to the classification given above, belongs to the second class. Moreover, the 
polarization obtained in this paper subtract the vacuum polarization at all A;, so that it vanishes 
identically in the vacuum. In contrast, the vacuum polarization of schemes A and B does not vanish 
away from the mass shell. For the 5 meson, it uses as in [ p4| ] the procedure of Kurasawa and Suzuki 
(KS) (thus belonging to the first class defined above) and then subtract the vacuum polarization 
at all k, Ilvac{M, k) = IIks{M, k) - IiKs{m, k). 

Several parameter sets were tried besides those given in the table of §4.1. In all cases, we find 
normal mode branches for the 5 and longitudinal p (see Fig. 1). For the /9, there are moreover two 
heavy meson branches. Contrary to the very similar a-uj system, no zero sound mode was found. 

Effective masses can be defined as usual as the solution(s) m* = w of the dispersion relation 
Z)(a;, A;) = at A; = 0. With the renormalization scheme A, the p meson mass increases with density 
and the 5 meson mass also slightly increases (see Fig. 2). With renormalization scheme B, the p 
meson mass first decreases, reaches a minimum at ~ 1.8 po ^^nd then slowly increases, whereas the 
delta meson mass increases reaches a maximum at ~ 0.8 po and then decreases. The p and 6 masses 
are almost constant with temperature (see Fig. 3) with both tti* and ms slightly decreasing. This 
figure was plotted using renormalization scheme A. In scheme C we reproduced the results of |44], 
with both the rho and delta masses decreasing as a function of density. 

This illustrates again the difficulties met with the standard renormalization techniques, which 
preclude a reliable prediction of the behavior of the effective meson masses in the medium. One 
cannot simply drop the contribution of the vacuum term by performing a normal ordering, since the 
contribution left out depends on the density through the effective mass of the nucleon. Moreover, 
pathologies appear in the dispersion relation (kinks, no clean normal modes) if one attempts to do 
so 1^ |5|, ^]. One could in principle apply a subtraction procedure at the one-loop order, however 
the various schemes used in the litterature ||2|, 38, 42| lead to widely different results. This is 
related to the fact that we are dealing with an effective theory which should enforce the scalings 
and symmetries of the underlying more fundamental theory, whereas the approximation made to 
the full many-body theory (here, RPA) blurs these concepts. There is some hope that one could 
solve the problem by applying "naturalness" and symmetry arguments [60, 
writing however, we must consider this a still unsolved problem. 



H]. At the time of 



One sometimes defines a mixing angle [62[ by 



1 

— arctan 

2 



2,/\^\W^ 



ml 



vat, 



(74) 



which is obtained from the diagonalization of the mass matrix for the mixed dispersion relation of 
two mesons A and B in the timelike region g^ > 0, 77^ < 0: 



M 




m 



B 



V'l^AB 



(75) 



det [q'^I - M] 



{q' -m\ + HaW -m%+ Ub) + v'^ab 
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The mixing angle is calculated at the solutions of the dispersion relation, there is therefore one 
for each branch. The mixing angle is represented for renormalization schemes A and B in Fig. 4 
as a function of density at vanishing temperature for a momentum k = 300 MeV. Schemes B and 
C yield values of 9spL of the order of a few degrees. The higher values found in scheme A results 
in fact no so much from the strength of the mixing, but rather from the fact that the difference of 
the effective masses of the 5 and p happens to be smaller in this scheme. 

The behavior with k is non monotonous. It is represented in Fig. 5 in renormalization scheme 
B at T = and ub = 3 rigat- We also calculated the effect of finite temperature. The mixing angle 
decreases with T in all renormalization schemes. This might reduce the effectiveness of S-p mixing 
as a mechanism invoked by the authors of [44| for dilepton production. 

4.3 Transverse ai-p mode and longitudinal ai 

Here we investigate the dispersion relation of the ai meson, and in what measure does the mixing 
with the ai meson modify the transverse p dispersion relation. 

The dispersion relation for the longitudinal ai mode is represented on Fig. 6 with renormaliza- 
tion schemes A and B. Only normal branches appear. 

The ai-p mixing vanishes at k=0. It therefore does not affect the effective masses, and m* 
coincides as calculated from the transverse or longitudinal dispersion relations. The ai effective 
mass is plotted on Fig. 7. With renormalization scheme A, it decreases slightly with increasing 
density. When using renormalization scheme B, a non monotonous behavior is obtained for the oi 
mass, which first presents a steep increase at low density, and then decreases again, in contradiction 
with the behavior expected from QCD sum rules [p3, p3, B5|. In the case of the ai meson therefore, 
the renormalization scheme A seems in better agreement with the results expected from other 
QCD-based models, whereas we saw that for the p meson, scheme B would have seemed preferable. 
A third renormalization scheme was also tried (let us call it scheme C), using the same procedure 
as for the 6 in [Q], (that is, subtracting the vacuum at all k from the expression obtained from 
the procedure of Kurasawa and Suzuki, see previous section and Appendix A. 2). In this scheme, 
the effective mass of the ai decreases. The behavior of m* as a function of temperature was also 
investigated; for all renormalization schemes it is almost constant as a function of temperature. 

At finite k, the ai-p mixing sets on. It does not appreciably affect the position of the normal 
branches. Nevertheless, it amplifies somewhat the (spurious) zero sound branch which may appear 
in the transverse part of the dispersion relation at high density. This mode corresponds to the 
p meson and is due to the high value of the tensor coupling fp/gp = 6.1 of the Bonn potential. 
It is weaker and appears at a higher density with the vector dominance value fp/gp = 3.7 and 
disappears completely if fp=0. Note however that a high value of fp is also supported by QCD 
sum rule calculations |£6|. This mode is present in both schemes a high enough density, but is 
stronger in scheme A. Such a mode would lead to divergences, or at least an enhancement of 
Friedel oscillations of the p contribution to the screened nucleon-nucleon potential [^ . It could be 
eliminated if a strong cutoff is applied to the p meson, of the order of 1200 - 1300 MeV. An other 



possibility is to use the contact term introduced in section §3.2, which we will need anyway for the 
pion. It was checked that the spurious zero sound branch appearing in the p dispersion relation is 
removed at all densities by the contact term used in next section. 

The strength of the ai-p mixing can be estimated by calculating the mixing angle, as explained 
in the previous section. In Fig. 8 it is shown at T = and k= 300 MeV. 

4.4 Pion mode 

If the dispersion relation is calculated with pure pseudovector coupling, and neither mixing with the 
oi meson nor Ansatze for short range corrections (a la Horowitz or Landau-Migdal contact term) 
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are taken into account, a zero sound mode is found already at saturation density in renormalization 
scheme A. The short range corrections must be added since such a mode is not observed. Before 
we pass to examine this point, let us first make a few more observations. 

Contrary to the non relativistic case, the zero sound mode is found to disappear again at higher 
density. Such a behavior was already noticed by Dawson and Piekarewicz [q^. The zero-sound 
problem is in part due to the high value of the cutoff A^ of the Bonn potential model, whereas 
several studies favour a lower value A^^ = 800 MeV. If we choose g'^/{A7r) and A^ = 800 MeV, there 
is only a tiny zero sound branch around k = 250 MeV, between ub = 0.72 ngat and ub = 1-2 nsat- 

In renormalization scheme B, no zero sound mode is found when pure pseudovector coupling 
of the pion is used without mixing with the ai. However this renormalization scheme appears less 
favorable when we consider the behavior of the effective pion mass. Whereas the pion mass stays 
approximately constant with scheme A, it strongly decreases with density when scheme B is used. 

In both cases, it is possible to adjust the effective mass at saturation density as measured in 
recent experiments ||6^ by adding a small admixture of pseudoscalar coupling { C 3 ip {iagT^^^ir .t 
+(1 — a)(/7r/wi7r)757^9;x7ff) V' + ^5o-7r7r"^o-c7r^ | ]69[ . Such an admixture has been considered from 
the theoretical |7^ as well as phenomenological |£|, ^ point of views. Gross et al. |7^ could allow 
a value as high as 25 % of PS admixture whereas Goudsmit et al. [^] found an optimal value of 
3.5 %. 

The crvr^ term must be introduced in order to compensate the contribution of the pseudoscalar 
coupling to the pion-nucleon scattering. The n — N scattering length can be reproduced when 
choosing ga-Kn = oP [g"^ / ga)^^^ . With this value, the pion mass can be adjusted to the value 
extracted from measurements of the pionic atom |^^ "T'7r(?^sat) = 1-1 i 0.03, with a = 0.0735 in 
renormalization scheme A and a = 0.118 in renormalization scheme B. 

As was to be expected, the zero sound branches are amplified when the mixing with the oi 
meson is taken into account, since it modifies the pion dispersion relation formally as would a 
Landau-Migdal contact term but with the opposite sign. Then it is present in renormalization 
scheme B as well. The zero sound branch can be eliminated as before by adding a contact term. 
The relevant formula when a PS admixture and a civr^ term are also present is 



n. 



ml - g^^^m^a + 11,^ = 



^) Upy + glU,s + 2g^^U 



m^ / m-n 



mix 



4 + '29a] (g^Kl^q, + ^KVq,' ' 



m^ I \ m 



g2+(4-25A|n 



ml 



PV 



UZ''q, = Il,y, nf/g, = ^np, (76) 

(For clarity we draw the coupling constants out of the definition of the polarizations in this equa- 
tion.) Mixing with the ai can somewhat increase the effective pion mass. This effect is less 
important after including the pseudoscalar admixture. 

The Horowitz Ansatz was also studied. The part of the dispersion relation which corresponds 
to the pion mode has the required form of a pseudovector coupling with a short range correction 
of the Landau-Migdal contact type. This removes the zero sound branch which appeared there for 
a pure pseudovector coupling. However, this Ansatz does also affect the rho meson mode through 
mixing with the fictitious "pion" introduced there. It is seen in dispersion relation Eq. (|4^ ) that 
there are also factors 1 + g'Tl-j^T, ^ + g'^TTL ■ It was found that these terms give spurious branches in 
the spacelike region whatever renormalization scheme was used. The most favorable case was that 
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of scheme A, which removes the branch due to 1 + g'Tl-^L entirely and part of that due 1 + g'Tl-KT 
at low density. These branches are an artefact of the unphysical components of the auxiliary 
pseudovectorial "pion" and should be removed carefully by hand at the end of the calculation. We 



therefore prefer using the contact Lagrangian of section §3.2, since this is free from such problems 
from the onset, moreover it permits easily to allow for a PS admixture, whereas Horowitz Ansatz 
can only be implemented for a pure pseudovector coupling. 

5 Conclusions 

This work gathers and extends results on the dispersion retations of mesons in relativistic nuclear 
matter at high density and finite temperature, as obtained in the RPA approximation to the 
quantum hadro dynamics. Besides the "standard" mesons a, oj, vr, p, results for additional mesons 
6, ai are given. Both mix with the p meson, the 6 with the longitudinal mode and the ai with 
the transversal mode. They therefore may represent an additional source of modification of the 
dilepton production at finite density and temperature. Also, new interest arises in the 6 meson in 
the context of the description of asymmetric matter in density dependent mean field theories |4G]. 



Various ways of introducing short range corrections in order to eliminate unobserved zero sound 
modes at saturation density were examined. It is seen that the ai meson mixes with the pion so that 
it acts exactly as would a Landau-Migdal contact term, but with the "wrong" sign. The Ansatz 
by Horowitz comes with the right sign, however it introduces spurious branches in the transverse 
channel which badly affect the dispersion relation in the spacelike region. A simple standard contact 
interaction in the Lagrangian does the job best. 

It was shown that the results are affected by the way the renormalization is performed in 
order to regulate the high momenta divergences. Without any renormalization, there appear kink 
structures in the effective meson masses as a function of density, no matter what a strong cutoff is 
applied. The expected branches are recovered when applying a renormalization procedure. Among 
the possible subtraction procedures, two schemes A and B were introduced and their predictions 
compared. Such important results as the behavior of the effective meson masses or the presence of 
zero sound modes differ widely whether one or the other scheme is used. In order to be consistent, 
a same scheme should be used for all mesons. Schemes of the A class have been used in previous 
litterature to study a, uj and it mesons. Scheme A yields results more in agreement with other 
models for the behavior of the vr, a and oi masses, but would predict an increasing p meson mass. 
On the other hand, schemes of the B class which give a decreasing rho meson mass, as favored by 
theoretical models and experimental data, would predict strongly increasing a and ai masses, and 
a rapidly dropping pion mass. This last result stands in strong disagreement with the expected 
behavior of a Goldstone boson and the recent experimental determination of 



A renormalization procedure which would respect the scalings and symmetries [30| of the un- 
derlying more fundamental theory of which the effective Lagrangian means to be a low energy 
approximation is clearly needed. 
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Appendix A: Polarizations 

There are contributions from the meson self interaction terms and from the particle hole insertions 
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The (retarded) particle hole insertions are given by 




n^^(p)= f d^pTr 


'r^s{p-^)T^s{p + ^) 


j^{p + i)-^ip-i) 

\ 2p.q - i e sign(po) 


where 







(77) 



S{p) 



-f.P + M, 



d 



{2ny 



■ 5(P2 _ m2) [9{poMp) + e{-po)n{p) - e{-po)] , 



d = 2 is the isospin degeneracy, and the F^ are the vertices, respectively for the a, uo, vr, p, ai 
couplings to the nucleon 



r^ G (5a,5L^7^,^57r75,^(/7r/m^)757^9M'5p7'' + wjj(^'"'Qiy,9aiJ5l'^) 



2M 



(78) 



As an example, we have 



Tr[{ga,i5r}GN{p-^){ga,l5Y}GNip+^)] 



4^[(V-M2 + ^)5M^ 



1 



k^'k'' + 2p^'p''] 



(79) 



Tr[{g,j^ - I^a^-k^}GN{p-^){9a,l,l''}GN{p+'l)] 



k. 



jp ^I'Oi, ^ ^ ' k. 



Tr[{ga,j5r}GN{p - -){gpY + j^a^'^k^jG^ip + ^ . 
^igaA9p + fp)e'""'^kaP(3 
In symmetric nuclear matter, the polarization matrix have the components 



(80) 



n,,. 












































^'-'qri 























^55 


^s, 

















n?. 


Wp 


^^ap 








































n^^ 



(81) 



In asymmetric nuclear matter the mixing polarizations II^p, H^p, IIo-5, H^i^, 11^^^ ... would not 
vanish. Asymmetric matter is discussed in a companion paper [pq]. 

The polarizations involving the ai meson are related to others by the following equations 



TT/iv ila Y\pv J ^a y TT-P^ ^ TT 



9l 



1 



Uir/m^, 



PV 

\2 rt'i "'"^'^ 



n 






^-J^Kn^^ 



These relations are used to simplify the expression of the pion propagator appearing in Eq. 



(82) 



n^ = n. 



ml -Ha 






ml /2 



(83) 
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A.l - Polarizations, real matter part 



The matter part of the polarizations is given by 



„2 
TT(mat) _ _ J yg 



2Xo+r--MV]T'o 



TT(mat) 



-d^Mq^ 






27r2 



n 



n 



(mat) 
LuL 

(mat) 



-d 



n 



(mat) 



-d 



27r2 
47r2 

A 

27r2 



7-2 _^o 



M^g^jO ^ (^2 _^ ^2) /2'2 ^ 2-o\ _ 4^^ji 



2x0+ ^-mV 4 



|-(mat) 
^5p 



n 



p 

(mat) 






nr"^^^ = -d ^^ ( 2M5p + ^g2 j ^2 



1^, 



T^ — —T^ 



2„2 



pL 



-d 



9pq 

27r2 



T2 -T^ 



"-i-n 



' 2m 27r2 






n 



(mat) 
PT 



, gp 

47r2 



M2g2j0 + (^2 ^ ^2) /j2 ^ jO\ _ 4^^jl 



f \ 2 „4 



2m/ 47r2 



M2j0+j2 jO 



^0 "T-^O 



<^(^)..^is 



2my^''27r2 



n 



n 



n 



(mat) 
vrPV 

(mat) 
a L 

(mat) 
a T 



-^f^)^ 



m^ 



-UaLL''''-UaTT^''-Ua 



fiu . 



2^2 



• 9aQ_ 

' 27r2 

47r2 



0° - (w^ + fc^) (Xo^ + if) + ^^kll 



M'^q^l[ 



n 



(mat) 
a Q 



'■'-ap 



d^^q'M^l', 

jd ^ ( „ _L 2 A/f ^ ^ fM^^a/?, 



'«??/3 






The integrals X^ which appear in these expressions are given by 



lm{p,iO,k) 



n—1 



dpp^Jp^ + M^] lra{p,uj,k) n(p) + (-l)™+"n(p) 



+1 



du 



(puY 



(84) 
(85) 

(86) 
(87) 
(88) 
(89) 



(90) 

(91) 
(92) 

(93) 
(94) 
(95) 
(96) 



{w^f{p^ + M2) - pku - q^/2){wy^{p'^ + M2) - pku + ^2/2) 



(97) 
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A. 2 - Polarizations, vacuum part 



The vacuum part of the polarizations comes from the —9{—p^) term in the definition of the rel- 
ativistic Hartree approximation to the Wigner function. It is divergent but can be renormahzed 
by dimensional regularization and subtraction of appropriate counterterms. The standard proce- 
dure is exposed in [§, ^]. In the case of the p and vr mesons, the derivative coupling makes the 
model nonrenormalizable in the usual sense. In practice however, a renormalization scheme can 
be defined at a given order of perturbation / loop expansion. One has the choice between several 
renormalization schemes, which were divided in ||5| into two subclasses: the "increasing rho mass" 
and the "decreasing rho mass" classes, in reference to the much debated issue of the interpretation 
of dilepton production in heavy ion collisions and the Brown-Rho scaling conjecture. To the former 
class belong earlier procedures by Chin or Kurasawa and Suzuki [37, 38]. In the latter class we find 



more recent procedures used by Shiomi and Hatsuda |^, Q. Both classes have advantages and 
drawbacks. In any case, a same scheme should be used for all mesons. 

To the Lagrangian (|^) we will have to add the following counterterm Lagrangian in order to 
perform the dimensional regularization and renormalization of diverging vacuum terms 



-^CT 



Z^{d^af + A^a^ + B^a^ + C^a^ + Z^Ff^'F^, + 



+ZaA>"'Apu +(^A, + Baa + Ca a^) a^ap +(F + G^a + H,cj' 



(a^d^TT) 



(98) 



Renormalization scheme A 

The first renormalization scheme is obtained by imposing that the polarization vanish in vacuum 
(M = m) on the mass shell g^ = m^, for each meson a = a, uj, 5, p, vr, ai, rj. We also require that 
the derivative with respect to q'^ vanishes on the mass shell for the mesons p and ai. For all mesons 



except the to, divergences of the type 



M' 



appear, which require for their cancellation 



counterterms like A^^a'^ + Ba-fJ^ + C^a^. The additional constants are fixed by imposing that the 
derivatives with respect to the a field vanish. 



dU 







d'^u 



shell 



da^ 







(99) 



shell 



The expressions given here are that of reference 

d gl 



for the a and to. We have: 



■Q(vacA) 



2 2^2 



6M2 log M - g2(log M + e) + AM'^O - (4 - ml) 6^ 
+(g2 - ml) {9„ - (4 - ml) 9,^) 
+(1 - M) (6 - m2 + (4 - ml)e„r. + 8 6,) 

Um + 80^ + 16( 



|(l-M)2(l8 + m2+(4 



m„ 



(100) 



where we express all quantities in units of the nucleon mass ("m=l") and with the definitions 

dx "2 



9{q^,M' 



y 



c^ + y)\/ x^ + 1 



with y = 1 



4M2 
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de 



de 



dM 



dH 



Jaq 



d{q^ 



q'^=Tn'^,M=m 



J- (vac A) . 



q'^=Tn^,M=m 



J (vac A) 



9M2 



q'^=m^,M=m 



(101) 



The expression for 11^ is identical to that for 11^ when replacing g^^, rritj, 9^ ••• hy gs, ms, 

9s ■■■ For the iv meson the renormalized vacuum polarization is given by 



TT/^i/(vac A) 



2 37r2 



2 r 



logM + ^ + 



2M^ 



1) 



1 



mf: 



[2 + ml) 9^ - 2 



q'g'"' 



fq''\w2) 



For the p we take the "scheme 2" expression of reference 

2M2 



TT/ii^(vac A) 
PP 



9l 



d 

2 1 3^2 



7r2 V 27n 



V2my 67r2 



log M + + ^(^ - 1) - ^ ((2 + ml) 9p - 2) 
M (logM + 9)-9p-{l- M) {ep + l + 9pm] 
6 M^ log M + 8 M^e + g2(iog M + 0) - (8 + m^) 0^ 



pJ^'p 



+(1 - M) (e + m^ + 16 ^p + (8 + ml) 9pm 



(1 - M)2 (l8 - ml + 16 9p + 32 0^^ + (8 + ml)9pmm 



-{q^-ml)[{8 + ml)9p, + 9p) 



q2gPU _ qf^qU 



(103) 



For the pion with pseudovector coupling, we obtain 

d 2 / /^ ^ 2 



(vac A) 



n 



T-r^(vac A) 



2 7r2 V"^7r 



M2(logM + 9)-9^ + il- M) (1 + 29^ + 9, 



^il-M)^{3 + 29^ + i9^m + ' 



iM 



^^^J^_jj(vacA) 
g2 (/tt/wI^) 



TTTT 



(104) 
(105) 



Since the polarization of the ai meson can be written as U^'^ = {ga/gu))'^Ilf^'^+ga{iT^TT/fn)'^^n^g^'^/q'^, 
it is tempting to apply the renormalization procedure with this decomposition. In this case, we will 
have Ilaa ^'^ = Hi {q^g^'^ — q^q^) + 112 9^^ where the expression of the Hi contribution is obtained 
from that of IliJ'^'^ by replacing g^, rn,^, 9^ by ga^ m^a, 9a- and that of 1X2 is proportional to IlTrT"^ 
when replacing f^jm-,^^ m-,^^ 9-,^ ... hy ga, m,a, 9a ... However, the 1X2 term is not orthogonal to 
q^q'^, whereas the polarization of the ai enters in the dispersion relations with the decomposition 
UaxT'^'^ + IlaLL^'^ + Ii.aQ q^q'^ /q^. It is therefore preferable to apply the renormalization procedure 

to n^a^™^) = li^aLT^a^" " Q^'f/'f) + ^^aQ^l^ff/o^- The 0^™^^ then happens to have the same 
structure as the polarization of the a meson. We obtain 

d 9l 






2 37r2 



GM^ log M - q^{\og M + 9) + AM'^9 - (4 - ml) 9a 
[q^-ml)(9a-{A-ml)9aq) 
+(1 - M) (6 - m^ + (4 - ml)9am + 8 0« 



+ 
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d2gi 
2 7r2 



-(1 - Mf (18 + ml + {A- mDOam^ + 80^ + 16^^ 

M2(log M + e)-ea + {l- M) (1 + 26a + 0am) 



AlU 



^fi^iy 



1 



il-Mfi3 + 29a + i9am + ' 



q^'q" 



(106) 



The vacuum contribution to the mixed polarizations 11^^, 11^ and 11^^ vanishes. 

Renormalization scheme B 

In this second renormahzation scheme, we notice that there appears in the polarizations di- 
verging contributions of the type M/e, M^/e, e ^ 0. Instead of cancelling this divergences with 
counterterms of the form a + ba + ca'^ , the original structure is preserved by subtracting only coun- 
terterms in the combination a^m — gaO') or a^m — go-c)^- The unknown constant o is determined 
by imposing n|sheii = 0. No conditions are imposed on the derivatives. 

Originally this scheme was motivated by the fact that it yields a decreasing effective p mass ^ . 
When applying it to the other mesons, we obtain 



-[-j-(vac B) 



d gl \ 
2 27r2 



GM^ log M + AM'^e - g2(log M + ^) - (4 - ml)e^ 



+ (M2 - 1) ((4 - ml)mle„q - 40,) + {q^ - ml) (o^ - (4 - ml)6„g)] (107) 



The polarization of the 5 meson is identical when relacing g^, m„, 0^ ... hy gs, ms, 6s ■■■ 

Since no divergences proportional to M or M^ occur in the polarization of the u meson, its 
expression is left unchanged by this renormalization scheme and its expression still given by Eq. 
(|102|). We recall the expression obtained for the p meson in scheme B (= scheme 3 of ||5|) 



TT(vacB) 
P 



d \g^ 
2 1 3^2 



lnM + 6' + 



2M2 



1 



,2m/ 67r2 



+ 7r2 \2m)^P 



q^ nip 

r'^^„ T\/r i c (I/tS/j , r,2 



A ((2 
m.'t \ 



+ ra'p)6p 



6M2 In M + 8M26I + k\\n M + 9) - {8 + m^) Op 
-{q^ - ml) ( (8 + ml) 9p, + 6p) + (M^ - 1) (m2(8 + m^) 9^ 



M{lnM + 9-9p) 



I {q'g^^'^ - q'^q''} 



For the pion with pseudovectorial coupling, we obtain 



-pf(vacB) 



d2^fU 
2 7r2 



m„ 



q^M"^ [\ogM + 9-9^] 



89, 



(108) 



(109) 



Finally, the polarization of the oi meson is given in scheme B by 

n..(vacB) ^ _nl7/)(5' 

d 9l 



nv 



q2 ) ^^aQ ^2 



n 



(vac B) 
aLT 



2 37r2 



n 



(vac B) 
aQ 



6M2 log M + AM'^9 - q^{\og M + 9) - {A - ml)9a 
+ (M2 - 1) ((4 - ml)ml9a, - 40.) + {q^ - ml) (6a - (4 - m^)^ 

2^4^'[logM + 0-0,] 
2 vr^ 



aq 



(110) 

(111) 
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A. 3 - Polarizations, imaginary part 

All imaginary parts can be expressed in terms of three integrals Ei, E2, E^. At finite temperature, 
for the calculation of the retarded polarizations, we have 
— for spacelike momentum: 



En 



dy 



y + 



0{y - yL)\n{y) - n{y + uj)j 

+(-1)" [y - 2 j ^(y - yu) My) - ^y - ^)} 



(112) 



— for timelike momentum with q^ < 4M^, the imaginary parts vanish: En = 0, 

— for timelike momentum with q^ > 4M^: 



En 



M 



UJ 



n— 1 



(-1)"'^ / dy[y--\ [e{y - y^) - e{y - yu)] {n(^ - y) + n{y) - 1} (113) 



with 



n(y) 



VL 



oPiy-f^) 



+ 1 



-1 



, n{y) 



ol^(y+t^) I I 



-1 



k\^-Lj 



, yu 



k^fK + uj 



A = 1-4 



M' 



In the following equations d is the degeneracy parameter [d = 2 for symmetric nuclear matter). 

^2\ 



Tm IIct = 

Zm I\-u)T = d 

Im U„^ = 

Im UpL = 

Im UpT = d 



.dA. ( m2 - ^ I El 
2TTk\ 4 / 



5^ 



27rk 
2 



9o. 



A-Kk 



-jEi - j^Es 



M^ + ^] i^l + |Ji?3 



i9cr9u} 2 Tit 771 



1 



2'Kk 



9l 



-jEi - j^Es 



IP „2j 



+ ..-^.^Mi., + ^ 



M2g2 El + I2S3 



1 



ATTk I ^f 



M' +'^)ei + |Ji?3 



fp ^2; 

''2m 



+ 2gpi^ q^M Ei 



2m J 



^mV + '{]ei-'^,e, 



ImUl" 



Imlll' 

Im liaL 



2TTk \mT^ J 
2Trk 



'*^'-t)^'4^= 
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Im YiaT 

Im liaQ 
Zm ^a-K 
Im liap 



Im J^f 



A-nk 



'^^^-t)^'4^» 



~27rk 



-d^^M^ E- 



d- 



9a I J TV 



a 

M 



+ l-a]M^ El 



2'nk Vmvr 
^ gg {gp + 2M{fp/2m)) ^ 

a : — 7-^ H -C/2 



47rA;3 



-l^-' 



Appendix B: Polarizations modified by contact terms 

In meson exchange models of the NN interaction, the potential is obtained by calculating the 
transition amplitude M = J2a,b=a,u;,n,p,s,ai,r,(U3'^aUi) G"^ (JhTaUi), with V" and G""^ beeing the 
relevant vertices and meson propagators. After taking the semiclassical limit and multiplying by the 
minimal relativity factors, one finally performs a Fourier transformation to obtain the potential in 
coordinate space. Besides the desired Yukawa-like terms with ranges characterized by the inverses of 
the meson masses, one also obtains a singular contribution at the origin 5{r) due to the assumption 
that all particles are pointlike. 

In a a-uj-Ti-p model, the singular contribution is given by |55] 



V{r) B 



al ^ al 



Am? 

,2 



2m' 



Im,'^ 



5{r) 



az 



+ - — -T1.T2 -\ ^ — K^ri.T2 



2m? Am? 



2m? 



^{^ ^<?l-'?2 



(114) 



It is not always enough to smooth this singularity by convoluting the result with form factors, 
especially in the case of the pion. It is then desirable to remove the singular piece by subtracting 
one or more contact terms. If this method is chosen, the dispersion relations of the mesons will be 
modified as follow. 

TT-ai-p sector 



When a contact interaction is introduced as 

the dispersion relations of the vr-ai^p sector are given by (^). The polarizations which appear in 
these expressions are given explicitely in this Appendix. 



II pL 
flpT 
f^aL 
flaT 



(ffp + 25ij^n^v)(5pn^v+&n^: 



2m (^^"-- + 2m "^0 ^ (1 - 2gnU^y) 

9, + 2,, An^,) Rpr + i2g,I^U-,)Ppr + ^ (,^n^, + An-, 



1 - 2gAU^^ 



^9aU^rA + (ap + 25iJ^n^,) RaT + {2gA^Ill^)PaT 



(115) 

(116) 
(117) 
(118) 
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PaT 

RpT 

RaT 
V 



(l - 2gRli 



(119) 



i;] ( gpiil^ + ^^l)j + 29flnl^ Ln^v + ^n^T 



1 - 2gRl[ly) gjilj, + 2gRgal(iyIily /V 



(l - 2gAliL) [gplily + ^n^x) + '^9Alil, [gplily + ^ri 
(i - 2gAnL) 5an^4 + 29Affan^^nL] /2? 

(l - 2gAnL) (l - 25^n^^) - AgAgRlllyli 



IV 



IV 



"-VA 



(120) 



with the definitions 11^^ 



..*n- + *(^)n: 






ll 

2m 



glUZ and H'^,'^ 



ap 



pv 

AT- 



al-uj sector 

Contact terms of the form 5{r) also appear in the contributions of the a and oj mesons exchange 
to the NN potential. Although the problem is less acute than in the case of the spin dependent 



part of the interaction mediated by the pion, one could wish for consistency [73| to remove these 
5{r) contributions by the introduction of further contact terms 

c 3 -gv{i^-fpip).{i^'j^i^) - gs{ipi^)-{i^il^) 

These terms modify the dispersion relations as follows 
[-q^ + ml]ai = glYlaOi + g^gu.^'^^ujip 



[-q^q"" + {q^ - rnDg^^i, = g.g^K^a, + gl (-U^tT'^'' - U^lL^") u^ 



with Ho- 









(1 - 2gvU^L)U„ - 2gvUl^r]^ 
(1 - 2 (71/ U^l) {l + 2gs n,) - Agsgv HL V^ 

gaQuj^'i^ 






(1 - 25y Yi^L) (1 + 255 n,) - ^gs9v nL ??' 



l-2gvTiu.T 



(1 - 2 gy n^ l) (1 + 2 55 n,) - 4 gs5y nL ry2 
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Fig. 2 Effective masses of the 6 and p mesons as 
a function of density, as calculated with renormal- 
ization scheme A (left panel) or B (middle panel) 
or C (right panel), at vanishing temperature. 
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Fig. 3 Effective masses of tlie 6 and p as a func- 
tion of temperature at three times saturation den- 
sity in renormalization scheme A. 
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Fig. 4 Mixing angle beween the 6 and p in the longitudinal mode at T = and k = 300 Mev, 

represented as a function of density, and calculated with renormalization scheme A (left panel) 
or B (right panel). 
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Fig. 5 Mixing angle beween the 6 and p in the longitudinal mode represented as a function of 

momentum (left panel) at ns = 3 rigat and T = 0, and as a function of temperature (right panel) 

at ns = 3 risat and k = 300 MeV. Both figures were obtained with renormalization scheme B 
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Fig. 6 Dispersion relation of the oi longitudinal mode, as calculated with renormalization 
scheme A (left panel) or B (right panel). The curves are labelled by the value of the density in 
saturation units n^/risat- 
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Fig. 7 Effective mass of the ai meson, as calculated with renormalization schemes A, B or C, 
plotted as a function of density at T = (left panel), and as a function of temperature T at 
ns = 3 nsat (right panel) 
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Fig. 8 Mixing angle in the ai-p transverse mode at k=300 MeV and T = 0, as a function of 
density. On the left panel with renormalization scheme A; on the right, with scheme B. 
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Fig. 9 Dispersion relation for the pion with pure 
pseudovector coupling, with and without mixing 
with the ai. The curves are labelled by the value 
of the density in saturation units n^/nsat and 
were obtained with renornialization scheme A. 
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